Strong correlation effects in 2D Bose-Einstein condensed dipolar excitons 
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By doing quantum Monte Carlo ab initio simulations we show that dipolar excitons, which are now 
under experimental study, actually are strongly correlated systems. Strong correlations manifest in 
significant deviations of excitation spectra from the Bogoliubov one, large Bose condensate depletion, 
short-range order in the pair correlation function, and peak(s) in the structure factor. 
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Two-dimensional (2D) dipolar excitons (DEs) with 
spatially separated electrons (e) and holes (h) are ex- 
tremely interesting due to increased lifetime, which per- 
mits to achieve different quasi-equilibrium exciton phases 
predicted for the system, e.g., 2D DE superfluid in ex- 
tended systems [H H, 0, 0, [1| and traps [1] , Bose-Einstein 
condensate , crystal phase [1, [§] and supersolid ■ 
A number of interesting phenomena, such as Josephson 
effects [m , light backscattering and other interesting op- 
tical properties [l^ can also be observed in the system. 

DE systems can be created in coupled quantum wells 
(CQWs) [H or in a single quantum well (SQW) applying 
a large external electric field (see theory in At 
present, a large experimental progress has been achieved 
in the study of 2D DE collective properties in CQWs 

and sqw 0. 

The superfluidity and coherent properties of equilib- 
rium e-h systems (in a sense of "spatially separated 
semimetal") have been also theoretically studied [H, d. 
An important progress was achieved by studying 
an electron bilayer in hi gh m agnetic field with 1/2 + 1/2 
filling of Landau levels [2l|. It can be proved that the 
properties of the system can be presented as a BCS state 
of a spatially separated (cornposite fermion) e-h system 
at zero magnetic field [l|, 0, The predicted e-h su- 
perfluidity and Josephson-like effects in this system were 
observed later on experimentally [23l |. It is worth notic- 
ing that the 2D dipolar Bose systems under considera- 
tion have been recently realized in atomic systems with 
large dipolar moments {e.g., for Cr atoms) [24| and polar 
molecules [25|. 

The majority of theoretical models describe the exci- 
tonic Bose condensate as an ideal or weakly correlated 
gas. Unfortunately, these approaches have a very limited 
region of applicability. Indeed, at small densities the re- 
pulsive dipole-dipole potential can be described by one 
parameter a, the s-wave scattering length, and the prop- 
erties are expected to be universal, i.e., to be the same for 
all interaction potentials having the same value of scat- 
tering length and, in particular, to be the same as in a 
system of hard-disks of diameter a. In fact, the latter 
is known to be weakly correlated only in ultra-rarified 
systems Q- So, the model of weakly correlated excitons 
holds only in ultra-rarified gases which have extremely 



low critical temperature. For real experimental excitonic 
densities such models can provide only a qualitative de- 
scription. Thus, a more accurate model should be worked 
out and a more precise study should be done in order to 
describe 2D DEs in quantum wells (QWs). 

This paper is devoted to a detailed microscopic study 
of 2D DEs by means of the diffusion Monte Carlo (DMC) 
technique. We prove that excitons are in fact strongly 
correlated in all the main up-to-now experiments with 
CQW [H, [13, [3 ill which low-temperature collective 
effects in exciton luminescence have been observed. We 
have obtained the following results supporting this fact: 

(i) The dimensionless compressibility C = 
{m^ / 27rh'^) / X and the contribution of dipole-dipole 
collisions to the chemical potential C ~ ficiri /2Tr h'^n 
are much larger than unity Ci C ^ 1 (x and fic are 
the corresponding dimensional quantities, m is the 
exciton mass, and n is the total exciton density in g^x 
spin degrees; gex = 4 in GaAs). For weakly correlated 
excitons one has <^ 1. 

(ii) The Bose-condensate density tiq of excitons at 
T = is 2-4 times smaller than the total density n (on 
the contrary, in the weak-correlation regime one assumes 
no « n). 

(iii) There are prominent peaks in both the pair dis- 
tribution function and the structure factor providing an 
evidence of short-ran ge o rder. Moreover, in the struc- 
tures studied in Refs. |l6l . [TtI . at typical exciton den- 
sity, e.g., n — 5 ■ 10^° cm~^ two peaks in the structure 
factor and three bumps of the pair distribution function 
are clearly visible. 

(iv) The excitation spectrum strongly deviates from 
the weakly-correlated Bogoliubov prediction although, 
even in the crossover between exciton and e-h regimes, 
there is still no roton minimum as found in Refs. [la . ITtI . 

m- 

(v) In the superfluid phase, the superfluid excitonic 
density n; is close to total density n 26]. In weakly cor- 
related systems, the quantity ni is significantly smaller 
than n at the superfluid transition. 

(vi) The superfluid transition temperature in an infi- 
nite system Tc is only slightly smaller than the degener- 
ation temperature T^eg, but the quasi-condensation tem- 
perature is 2-2.5 times larger than T^eg [23]. According 
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structure 


D, nm 


1/4, 


-2 

cm 


Ref 


GaAs/AlGaAs CQWs 


15.5 


1.56 • 


10^° 


[14J 


InGaAs/GaAs CQWs 


10.5 


1.83 ■ 


10^^ 


[15] 


AlAs/GaAs CQWs 


13.6 


2.64 ■ 


IQio 


[16] 


GaAs/AlGaAs CQWs 


12.3 


3.94- 


10^° 


[17] 


GaAs/AlGaAs CQWs 


14.1 


2.28 ■ 


IQio 


[18] 


GaAs SQW 


^ 6 


^7. 


10" 


[20] 



TABLE I: Experimental structures and values of their e-h sep- 
arations D and dimensional densities n = I/sq correspond- 
ing to dimensionless density n = 1. The exciton mass is 
m = 0.14mo for InGaAs/GaAs CQWs and m = 0.22mo for 
other CQWs and SQW, with mo being the free electron mass). 



to the theory of weakly correlated excitons these temper- 
atures are logarithmically small compared to T^eg- 

(vii) The excitonic Bose condensate profile in a 2D 
large-size harmonic trap at T = differs appreciably 
from the Thomas-Fermi inverted parabola. 

For Timofeev's experiment with SQW we find ex- 
citons to be intermediately correlated. 

We have performed DMC simulations at T = of a 
2D homogeneous system with N = 100 DEs in a square 
simulation box with periodic boundary conditions and 



for 12 different dimensionless densities n = nxl 
where i is an integer number (1 < i < 12), and 



= 2' 



md 



Xq 



> 



(1) 



is the characteristic length scale for the excitonic dipole- 
dipole interaction V{r) = f^Xo/mr^. Here, d = eD > 
is the exciton dipole moment. 



D = 



{\^,{z)\'-\Mz)\')zdz 



is the effective e-h separation, with ipeiz) and i^hiz) the 
wavefunctions of e and h in the Oz direction, e > is 
the hole charge, e is the dielectric constant of the QW 
structure (e « 12.5 in GaAs), and eo ~ l/An is the per- 
mittivity of vacuum. 

Parameters of recent experiments [3 [1^, [IB, [13, [3, 
[20I I are summarized in Table [l] 

We note that the excitonic collective properties are ex- 
perimentally observed in Refs. [IB, [13, [3] at dimen- 
sionless densities l/4<n<4. At the lowest density 
n ~ 1/256 of our computation the su per fluid crossover 
temperature in the DE system of Ref. 20] is rather low, 
Tc « 0.12 K. At the highest density we have considered, 
n = 8, the dimensional density of DEs in the structure 
of Ref. [1^1 is rather large, n « 1.2 • 10^^ cm~^. 

The trial wave function used for importance sampling 
in DMC is similar to the one used in our previous stud- 
ies; for details of the calculation we refer the reader to 
Ref. [H. 

For each density, we compute the ground state energy 
per exciton Eq/N (see Fig. 1(a)). The DMC results are 
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FIG. 1: (a) Ground-state energy per exciton Eo/N as a func- 
tion of the density (we use h = m — xo ~ 1 units; length 
scale Xo is given by Eq. HI). Squares, DMC; solid line, Eq. 
0. (b) Dimensionless energy per exciton (solid line), colli- 
sional contribution to chemical potential (dashed line) and 
compressibility ( (dotted line) calculated from Eq. ((2]) as a 
function of the density. 



reproduced within a 0.025% accuracy by the polynomial 
fit (now we go to /i = 771 = = 1 units, so, n = n) 

Eq/N ~ fle exp(fee In77+Ce In^ n+d^ In'^ n+e^ In^ n), (2) 

with Qe = 9.218, be = 1.35999, Ce = 0.011225, 4 = 
-0.00036, and = -0.0000281. 

Using Eq. ^ we calculate the coUisional contribu- 
tion to the chemical potential fic — dEo/dN 29] and 
the inverse compressibility 1/x — d^c/dn. The re- 
sults for Eq/N, fic and 1/x in terms of the quantities 
C" = 2E0/NT0, C' = fic/To and C = n/xTo are shown in 
Fig. 1(b). (Here Tq = QexTdeg — 27r7i is the degeneration 
temperature of the spin-polarized DEs and Tdeg is the de- 
generation temperature of real (spin-depolarized) DEs.) 
We find that at densities 1/4 < 77 < 4 the dimensionless 
compressibility C is 3 < C ^ 8, i.e., C, ^ 1- This is an 
evidence of the presence of strong correlations between 
2D CQW DEs whose collective state was investigated in 

Refs. 0,[i3[ill- 

In Fig. 2(a), we show results for the polar- angle- 
averaged one-body density matrix 

/.27r 

pi(r) = / (*+(r)^'(0))d(^/2^. 



Here ^'(r) is the exciton field operator and the brackets 
(...) denote ground-state averaging. 

The asymptotic behavior of ^1(7") matches the result 
derived from a hydrodynamic calculation (see Fig. 2(d)) 
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FIG. 2: One-body density matrix (a) and pair distribution 
function (c) for the system of = 100 excitons at 12 den- 
sities (n = 2'~®, i = 1,2. ..12). (b) Bose-condensate fraction 
and polynomial fit Q as a function of the density, (d) One- 
body density matrix pi(r) shown in (a) at density n = 1 and 
large distances r: dash line, DMC simulation; thin solid line, 
hydrodynamic calculation. 
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giving us additional confidence on the accuracy of the 
DMC simulation. 

In Fig. 2(b) we plot DMC results of the Bose- 
condensate fraction no/n and the corresponding polyno- 
mial fit on top of the data, also in terms of Inn, 



no/n — flg exp(6Q Inrt 



'In^n + d'^ln^n), (3) 



with parameters ajf = 0.3822, = -0.2342, cjf = 
-0.02852, and d'^ = -0.001594. The results reveal a 
strong Bose condensate depletion at densities 1/4 < n < 
4. This is again a clear manifestation of strong correla- 
tions between 2D CQW DEs. 

In Fig. 2(c), we show the results of DMC simulations 
for polar-angle-averaged excitonic pair distribution func- 
tion 

P2{r) = r {p{r)m)dv/27r 



(/5(r) — ^'+(r)^(r) is the excitonic density operator). 
The prominent bump in the pair distribution function at 
densities 1/4 < n < 4 is an evidence of both the short- 
range order and the strong correlations present in DE 
systems. 

The structure factor S{p) is related to the Fourier 
transform of the pair distribution function P2{'r). We 
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FIG. 3: Exciton structure factor (a) and excitation spectrum 

(c) for 12 densities (n = 2'"^, i = 1,2. ..12). (b) Sound veloc- 
ity calculated from the slope of the structure factor (squares) 
and from Eq. ([2]) (solid line) as a function of the density. 

(d) Comparison of excitation spectra at different densities 
(n — 2'"', i = 1, 2. ..12) with the Bogoliubov prediction. 



plot the DMC results of Siji) in Fig. 3(a). A sharp peak 
in iS'(p) at densities 1/4 < n < 4 is consequence of the 
short-range order in 2D CQW DE systems. 

Using the DMC results of the structure factor 5(p) one 
can predict an upper bound for the excitation spectrum 
by using the Feynman formula [131 



25(p) 



(4) 



(see Fig. 3(c)). We see that the collective excitation spec- 
trum of the 2D DE system is rather far from the Bogoli- 
ubov form (see Fig. 3(d)) in the range 1/4 < n < 4. 
This fact suggests the relevant role of correlations in the 
excitonic collective state in CQWs. 

We compare two calculations of sound velocity Csi 
from the low-momentum structure factor slope (cg = 
p/2S{p —> 0); see Eq. (|4])) and from the compressibility 
1/x {cs = \AVx)j where x is calculated by differentia- 
tion of the energy fit ([2]). The good agreement between 
the two calculations (see Fig. 3(b)) is expected from an 
exact calculation like the present one and supposes a cru- 
cial internal consistency check. 

In Fig. 4(a), we show results derived from the Lan- 
dau formula of the tem pera ture dependence of the local 
(vortex-unrenormalized [3lL[3^.[33l|) superfluid fraction of 
2D quasi-condensed [IJl DEs ni (T) /n far from the quasi- 
condensation crossover. Results for the Berezinskii- 
Kosterlitz-Thouless (BKT) superfluid transition temper- 
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FIG. 4: (a) Temperature dependence of the local excitonic 
superfluid fraction in g^x = 4 spin degrees at 12 densities 
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2'- 



1,2. ..12; solids). Dotted lines depict the for- 



mal extrapolation of the calculation to the quasi-condensation 
crossover region, (b) Quasi-condensation temperature and (c) 
BKT transition temperature in the infinite system as a func- 
tion of the density. 
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FIG. 5: (a) Total and Bose condensate profiles of harmoni- 
cally trapped excitons in LDA at Ce = 10 for 10 densities in 
the trap center (n(0) — 2'~^, i = 3, 4. ..12). (b) Comparison 
of the Bose-condensate profiles at different n(0) with the TF 
inverted parabola. 



ature [H, [H, in a 2D infinite system Tc as well 
as the quasi-condensation temperature Tq are shown in 
Figs. 4(b,c). 

We see that at densities 1/4 < n < 4 the local super- 
fluid fraction ni/n at T < Tc is close to unity. Moreover, 
the quasi-condensation temperature at 1/4 < n < 4 
is 2-2.5 times higher than the degeneration temperature 
Tdeg = To/gex = 0.25To and the BKT superfluid tran- 
sition temperature Tc ~ 0.22ro is only slightly smaller 
than Tdeg- AH these features are manifestations of the 
strong correlations between 2D CQW DEs. Besides, our 
results show that in real, strongly correlated, DE sys- 
tems quasi- condensate and superfluid phases are much 
easier to achieve experimentally than if DEs were weakly 
correlated. (For weakly correlated DEs temperatures Tc 
and Tn are logarithmically small in comparison with T^eg 

Mm)- 

Based on our DMC simulations of homogeneous infi- 
nite systems (Eqs. ([2]) and Q) and applying local den- 
sity approximation (LDA) we have calculated the total 
n(R) and Bose-condensate rto(R) exciton profiles in a 
harmonic trap of large size {^/n{0)L'^p > 1) at T = 0. 

Here, Ltf ^ ^/2NJ^m{fi) is Thomas-Fermi (TF) radius 
and ^ 1 is number of DEs in the trap. The results 
obtained for both functions are shown in Fig. 5(a). 

If the electrostatic contribution is large = A^e/7o = 
2e'^D/e — 10 [11] (thi s is its typical value for the ex- 
periments of Refs. [IGI . [TtI [TsI). the total profile is in 



very good agreement with the TF inverted parabola (see 
Fig. 5(a)) [S^. However, the Bose condensate profile at 
1/4 ^ ^ 4 differs essentially from the inverted parabola 
(see Fig. 5(b)). The latter fact shows a failure of the 
theory of weakly correlated excitons for which the differ- 
ence between Bose condensate profile and TF inverted 
parabola is logarithmically small Q • 

Correlation effects in a dense 2D DE gas in CQWs can 
be detected by observation of the exciton luminescence in 
in-plane magnetic field. Indeed, in absence of magnetic 
field, the momentum of a recombined exciton p accord- 
ing to the momentum conservation in one-photon exci- 
ton recombination is equal to the QW plane projection 
[(Su/cq) sin 6] of the momentum of the emitted photon 
kuj/co [33] ■ Here is the angle between the normal and 
emitted photon in free space, cq is the light velocity in 
free space, and lo is the photon frequency. However, if 
there is in-plane magnetic field , the dispersion curve 
is shifted by the quantity pn = eDH^^/co [i^. This 
results in the following connection between the photon 
angle 9 and the exciton momentum p: 

(ftw/co) sine = |p - ph|. (5) 

If one considers only normal luminescence (6 = 0), then 
P ~ Ph (0), and thus Sp — Sp^. Hence, for spectral- 
angle luminescence along the normal in the field one 
obtains 

lf^,iu)^lf^,Siu;-n + epJh) (m»PT), (6) 
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where /|^q is the spectrally integrated angle lumines- pr ^ eDH\\/co ^ h^/n in ([7]) corresponds to the field 
cence along the normal in the field , Q. is the exciton 
resonance fi:equency for luminescence dX H\\ =6 = 0, 
and pt = T / Cg is the typical thermal momentum of Bose 
condensed DEs which is the boundary between thermal 
p ^Pt and zero-temperature pr ^ p <C h\/n regions of 
quasi-condensate phase fiuctuations [34| . 

From Eq. we find that the observation of the 

magnetic-field-induced normal excitonic luminescence 
line shift enables one to directly measure the excitation 
spectrum Sp. In the structure of Ref. [i^ {D = 12.3 nm) 
the density of the exciton gas is n = 2 • 10^° cm~^, the 
field range < iJii < 8 T corresponds to the momentum 



range {c^T / eD) m / (T^^ < H\\ < c^h^/eD, or, 0.1 T 
< H\\ < 0.75 T. 

The Bose condensate luminescence at T = entering 
intoEq. ^ can be approximately calculated as thermal- 
phasc-fluctuation quasi-condensate luminescence at low 
temperatures 



= K TLq 



- 

^0=0 



To 



(T « T,), (8) 



range < ph/^, < 1.5-10 cm^ . This range covers 
both long-wavc-lcngth (hydrodynamic) and intermediate 
scales up to short-wave-length (ideal-gas) ones on which 
^PH ~ Plf/2"T'- Note that the thermal momentum pT at 
T = 0.2To = 1 K < Tc,Tq, m = 0.22mo and C = 3.7 (see 
above) corresponds to a field H'^ = {cqT / eD) y/ m / (Tq « 
0.2 T. 

By measuring the slope of the excitation spectrum 
one can obtain sound velocity Cg ~ Spn/PH {ph ^ hy/n; 
see Fig. 3(c)). The dimensionless exciton compressibility 
is then found as 

Alternatively, one can determine the dimensionless com- 
pressibility C by measuring the magnetic field dependence 
of spectrally integrated angle luminescence along the nor- 
mal 



[pT < eDHu/co < hy/n). 



. (7) 

Here qr = huj/cQ ~ HQ/cq is the radiation zone in free 
space, Iq = K^hQ / TQ)no is the Bose-condensate lumines- 
cence at T = 0, To is the lifetime of a spin-depolarized 
exciton at T = in pure CQW structure, and k = 1/2 
(k = 1) if the luminescence is measured on one (both) 
sides of CQW plane. 

At n = 2 • 10^° cm-2, D = 12.3 nm, T = O.ITq = 
0.5 K, m = 0.22mo and ( = 3.7 the momentum range 



Finally, the measurement of the quantity ([8]) enables 
one to determine an important parameter of the exciton 
interaction, the Bose-condensate density at T = 

^o--^/o. (9) 
In conclusion, Bose-Einstein condensation, superfluid- 
ity, and microscopical properties of 2D dipolar excitons 
in SQW and CQWs have been studied by ab initio sim- 
ulations and analytical calculations. In all the recent ex- 
periments on excitonic Bose condensation in CQWs the 
excitons are strongly correlated. We have shown that in 
a strongly correlated regime the excitonic collective state 
is much easier to be achieved than in weakly correlated 
DEs. Our calculation of the microscopical properties of 
2D DEs enables one to investigate quantitatively QW 
excitonic Bose-Einstein condensation, superfluidity, lu- 
minescence and nonlinear optical effects by means of the 
hydrodynamic method in quantum field theory. 
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